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We consider the Casimir interaction between two spheres corresponding to massless Dirac fields 
with MIT-bag boundary conditions. Using operator approach, we derive the TGTG-formula for the 
Casimir interaction energy between the two spheres. A byproduct is the explicit formula for the 
translation matrix that relates the fermionic spherical waves in different coordinate systems. In the 
large separation limit, it is found that the order of the Casimir interaction energy is L~ 5 , where L 
is the separation between the centers of the spheres. This order is intermediate between that of two 
Dirichlet spheres (of order L~ 3 ) and two Neumann spheres (of order L~ 7 ). In the small separation 
limit, we derive analytically the asymptotic expansion of the Casimir interaction energy up to the 
next-to-leading order term. The leading term agrees with the proximity force approximation. The 
result for the next-to-leading order term is compared to the corresponding results for scalar fields 
and electromagnetic fields. 
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I. INTRODUCTION 

The Casimir effect induced by the vacuum energy of a quantum field is an important topic of research in various 
areas of theoretical physics such as quantum field theory, gravitation and cosmology, atomic physics, nanotechnology 
and condensed matters. In the end of the last century, the advancement in Casimir experiments have driven the 
acceleration of the theoretical research of Casimir effect. In the beginning of this century, we have observed a flourish 
in the Casimir research into new eras. In particular, there is a breakthrough in the understanding of how to compute 
the Casimir interaction between two nonplanar objects using multiple scattering or mode summation approach after 
the hardwork of a few groups of researchers [M3. 

In these ten years, many research have been done to understand the Casimir interaction between two nontrivial 
configurations such as a sphere and a plate, two spheres, two cylinders, a cylinder and a plate, etc. Using multiple 
scattering formalism, one can show that the Casimir interaction energy can be written as the integral over the trace 
of the logarithm of a matrix, which is the multiplication of four matrices, two of them are the T-operators of the two 
objects, and another two are the translation matrices that relate the coordinate systems used for the two objects. 
The two T-operators are intimately related to the scattering matrices of the objects and can be easily computed by 
matching boundary conditions. The difficult part of the problem is to derive the translation matrices, and one usually 
needs advanced mathematical tools. 

In the pioneering work [2(J (2l| , Milton has considered the fermionic Casimir effect of a spherical bag. Subsequently, 
Casimir effect of fermionic fields has aroused considerable interest f22l4ril . Nevertheless, to the best of our knowledge, 
no one has considered the fermionic Casimir interaction between two nontrivial objects. In this work, we take the 
first step in this direction. We consider the fermionic Casimir effect between two spheres exterior to each other and 
compare the results to that of scalar fields and electromagnetic fields. The fermionic field satisfies the Dirac equation 
and is subject to the MIT-bag boundary conditions. As has been well-known, although fermionic fields have spin 
and statistics that are different from that of bosonic fields, they also induced attractive Casimir force on two parallel 
plates. Hence, it is natural to expect the same result in the case of two spheres. We show that this is indeed the case 
by computing the large separation and small separation asymptotic behaviors. 
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II. THE CASIMIR ENERGY 

Consider two spheres A and B outside each other with radii ta and tb respectively. Let L be the distance between 
the centers of the spheres and let d be the distance between the spheres. Then d = L — rA — r’s- 

We want to study the Casimir effect due to the vacuum fluctuations of a massless spin 1/2- fermionic field with 
MIT-bag boundary conditions on the spheres. The equation of motion is the Dirac equation 

i^V ^ = 0, 

where = d M TT^, and T^ is the spin connection. On the boundaries of the spheres, MIT-bag boundary conditions 

are imposed: 


(1 + 


boundary 


= o, 


(i) 


where is the unit outward normal vector. 

To derive the formula for the Casimir interaction energy between the spheres, we use the formalism we developed 
in [19]. The TGTG formula for the Casimir interaction energy is given by 

E Cas = - — / In det (I - T A G AB T B G BA ) . (2) 

J o 

Here £ is the imaginary frequency. The minus sign in front of the integral appears because we are dealing with 
fermionic fields. In this formula, T A and T s are respectively the T-operators of the spheres A and B, <G AB and G BA 
are the translation matrices that relate the bases of wave functions in different coordinate systems. 

In spherical coordinates, the wave basis are parametrized by ( j,m ) with j = 1/2, 3/2, 5/2,... and m = —j , — j + 
1, • • ■ ,3 ~ b 3 - The y are given by 


C,i* - 


=c; e 


iujt 


«/,(=•=)>* iujt 
t jm,2 j e 


V =F if* +h (kr)n jj+h , 
' fJ+±( kr ) n jJ+bm 


where k = lo/c, and ip( ' are respectively the positive energy modes and negative energy modes, * = reg or out 
for regular or outgoing waves, 


trw =ji( z ) = 

fr\z) =hf\z) = 


Qjim are the spherical harmonic spinors fid : 


n 


jlm 


n 


jlm 


( hi + 


J + m 

2 j 


Y, 


l.m— Tv 


J ~ m Y 

VV 2 j r km+^J 


( _ j ~rn+ 1 \ 

V 2j + 2 l ' m - 


j + m+1 

V V 2 j + 2 l ’ m+ y 




Cj es = i- j+ i, C° ut 
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and the constants 
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are introduced to facilitate the change to imaginary frequencies. 

Matching the boundary conditions on the spheres, we find that the T-operator for sphere A is given by 


/ rp(±),A 

'tp(±)> j 4 _ / ^ jm, 1 

jm l 0 


T, 


0 

(±)M 

jm ,2 / 


with 

T (±),A = Ij{nr A )TiI J+1 (Kr A ) 

Kj{^r A ) ± iKj +1 (nr A ) ’ 
t (±),a = Ij(Kr A )±iIj+i{nr A ) 

Jm ’ 2 Kj(nr A ) =F iI<j + i{nr A ) ’ 

and k = in. The T-operator for sphere B T B is given by the same formula by changing r A to rg. 

The derivation of the & AB and G BA matrices are more complicated. Their components are defined by the relations 




( ,/.(±),°ut 
t jm,l 
/(±),out 
\'j rn , 2 


=YYg ab ’^\ 


f ,/>( ± )’ re g 

t jm,l 


(x,w) 


^•Sr s (x,w) y 


(x,w) 


-Wr BA,(:l: ) I L S ( X ') W ) 


(3) 


where x' = x — L, L = (0,0, L). 

To compute these matrices explicitly, let us use the operator method introduced in |43} and developed in 0. As 
in 0,0 , define a differential operator Vim by 


Vim =(-l) 
Vl-m = 


121 + 1 (l — to)! / d x + id, 
4-7T (l + to)! \ ik 


P, 


(m) 


21 + 1 (l — to)! / d x — id j 


ik 


p(m) i - 

' 1 ik 


47t (l + m.y. 

Here \z) is the m times derivative of the Legendre polynomial Pi(z). This operator is such that 

Vi m e tkr = Y l m(0kAk)e^ r , 


where 


k = (fci, & 2 , ) = (k sin 9k cos (f>k , k sin 9k sin cj>k , fc cos 9k ), 


and 


r = ( x,y,z) = (r sin 9 cos </>, r sin 9 sin (j> , r cos 9). 

It has been shown in 0,0 that 

Pimjoikr) =i l ji{kr)Yim(9,4>), 
Vi m h^\kr) =i l hf\kr)Y lm (6, 0); 


and 


•1 n 27T /»7T 

jo(kr) =— / d(j>k / d9 k sin^e^' 1 
4 tt Jo do 


4 1} 0) =^: 


dk x 


f>oo ik x x-\-ik y y±iyJ k 2 — k 2 — k 2 z 


dky 


k k^ — k2 — k^ 

IX\ / IX lx x IX y 


z^O. 


(4) 


( 5 ) 
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From (|4]), we have 


Now, define 


CJi* =C*e^H~^ 


2 j 




13 - m 
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Using the orthogonality of the spherical harmonics and the formulas ©, one can check directly that 

oTiuit 


Vf'm', i-® reg (x,a;) 

^L',2 •v , S£ , r s ( x .w) 


=(-!)■ 


x=0 


47T 


~$jj' $mm' i 


= 0 , 


x=0 




(±),reg. 


=0 


x=0 




i e^ iut 


x—0 


_( 1) J+ 2 j ^ Sjj'Smm’ 


From these and the definition ©, we find that 


G A ' B ’ ( ; ±) , =(-iy~hTre ±iu,t 

im.i'm' V / 


^m.l ■ ^)| x=0 - iV fm,2 ' ^( X '’ W ) l*=0 


I ^31 J ,OUt / / \ 

W'rn' 2 ( x > w ) 


jm,l rj'm',2 


x—0 2 


I J jOUl / / \ 

Tj'm',2 ( X > W )| 


After some computations, we obtain the following formulas for and Gy A ’/^: 


G AB ’ ( , ±} , =(-iy'+ m (5 m „ 

jm,jm V / ^m,r, 


7T (j - - to')! /' c 

2 V 0 + m)!0' + m / )! i 0 


dO sinh 9e~ K,L coshe 


0 + m)P m _ i 2 (coshf?) P™^ 2 (cosh0)\ /(/ + m')P™_ i 2 (cosh#) i(j' — m' + 1)P™ i 2 (cosh#)^ 
2 3 2 II •? 2 J “r 2 


m+r 


*0 ~ TO + 1)-P + i 2 (cosh#) iP. + 


(cosh 9)1 \ —P” ! _^ 2 (coshfi 1 ) 

' ' 3 2 


iP^ + i 2 (cosh#) 


and 


r j'm' ,jm 


1 0 \ ~AB,(±) /I 0 
0 -1 I ^jrnj'm' l 0 —1 


Hence, after some simplification, the Casimir interaction energy between two fermionic spheres can be written as 


Peas — 


dflndet(I-M). 


( 6 ) 


where 



0 


rp-\-,A 

1 jm,l u \ V^rl 
n r +A Y im,r 

U -Urn,2/ y, 


T~ 


+,B 


j"m, 1 

0 


T. 


+ ,B 


V: 




j m,j mi 


with 


V 


' 0 ~ ra)!Q" ~ to ) ! f 

(j + m)\(j" + m)! Jo 

0 + m)P m _ i 2 (cosh#) 


d# sinh #e 


-kL cosh 0 


P™t 2 (cosh6»)' 

J 2 i 

, — 0 — m + 1)H!“ i 2 (cosh#) Pj"t 2 (cosh#) ; 


0" + m)P .„_i (cosh#) (j" — m + 1)P., |(cosh0) 
•? 2 •? 2 


-P”( A i(cosh0) 
J 2 


p 2+l(cosh0) 
J ' 2 


and 


r2 _^rl 


In ©, we have multiplied a factor of 2 taking into account the contribution from positive and negative energy modes. 
Although the components of the T-operators of the spheres T" 4 and T s are complex, one can check directly that after 












6 


taking the trace, the contribution to the Casimir interaction energy from positive energy modes and negative energy 
modes are both real and equal to each other. 

We would also like to remark that using the identities of associated Legendre functions, one can show that the 
product 


( (j + m)P m _ i 2 (cosh0) (cosli0)\ f (j" + m)P™_\ (cosh 6) (j" — m + 1)P™ ? (cosh#) 

3 2 i 3 2 i { 2 i 3 +2 

(j — m + 1)P™ i 2 (cosh0) P™i 2 (cosh0) / \ — P™^?(cosh0) P™^|(coshd) 

\ 3 i 2 3 ' 2 / \ 3 2 3 <2 

is a matrix of the form 

A B\ 

-B -A) ' 

In fact, V] -/ m can be expressed as linear combinations of Kjn (kL) with coefficients depending on and m. 

However, it would not help our computations later. 

In the following, we are going to explore the asymptotic behavior of the Casimir interaction energy when the 
separation between the spheres is large and when the separation is small. 


III. THE LARGE SEPARATION ASYMPTOTIC BEHAVIOR 

When the separation between the spheres is large, i.e., when L r A ,r B , the leading contribution to the Casimir 
interaction energy comes from terms with lowest j and m, namely with j = 1/2 and m = ±1/2. Replacing n by k/L, 
we find that the leading contribution is 


P^Cas 


2 he 
nL 


dn trMo 


(7) 


where 






The factor 2 in front of the integral in © comes from m = 1/2 and m = —1/2 which give equal contributions. When 
L » r Al r B , 


Ti 4 = T 2 a 



2 i / nr A \ 2 

7 VT) ’ 


and a similar expression for and . 


Hence, 


V£=V 


2 

0 


dd sinh de~ K cosh 6 


1 

— cosh 6 


coshd\ 

1 J 


= J_ ( K K + l\ K 
ft 2 ( K + 1) — K J 


Pc as 


2 hcr\r B 
7 tL 5 

Ahcr 2 A r 2 B 

nL 5 



6 hcr\r B 
nL 5 


/c±l\ 

~ K J 



-(« + l) 


e_2K 


Namely, in the large separation regime, the order of the Casimir interaction energy is L 5 
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For two Dirichlet spheres, the large distance asymptotic behavior of the Casimir interaction energy is 


E 


D 

Cas 


hcr A r B 
47 tL 3 ’ 


whereas for two Neumann spheres, 

N 161 hcr\r% 

Cas ~ 9 6nL 7 ’ 

and for two perfectly conducting spheres, 


e: 


Cas 


143 hcr\r% 
167rL 7 


Hence, we see that in the large separation regime, the order of the Casimir interaction is intermediate between the 
Dirichlet and the Neumann case. 


IV. THE SMALL SEPARATION ASYMPTOTIC BEHAVIOR 


The derivation of the small separation asymptotic behavior is more complicated. Based on the pioneering work 0 , 
a perturbative machinery has been developed in [ll|j 0^0. We use the same method with necessary modification to 
the present scenario. 

First we expand the logarithm in 0 and get 


he ^ 1 


£ Cas = -£ 

TT ^' 


dK E 


E 


E 


S—o 0 m=±!,±f,....7o = M,|m| + 5 >... j a =\m\,\m\+±,.. 




with 




TT 2 (T+’ A 


0 T, 


+ ,A 


jim ,2 / 


Vv 1 


T+' b o \ 

| i>,i u \ V 2 

fi m l n T + ’ B 7 


Making a change of variables 


£ = 


C4 , rB / . x 

a= -, b = -— r —, to = n(r A +r B ), 


ta + vb' r A + rs' ta + tb 

jo = 1, ji = l + m 

b b b 

Ji — z~ [h + h+ 1 ) + 91 — — l + z- \P-i + n-i+i) + qi , 
Za a Za 


uj = 


l\J\ — T 1 


and replacing summations by appropriate integrations, we have 

Eg 


J%C °° 1 /* 1 (Z'T" f °° f °° f 00 f 00 S 

^ . i i / — 5 / , 9 / dll / dm / dn\... / dn s II (8) 

r* ^ I -*■ J 0 T V 1 T J 0 J —OO — oo J — OO • n 


7T7\4 


s—0 


2=0 


Next we expand each term in small e keeping in mind that l ~ e 1 , n*, ~ e 2 , r ~ Using the formula 

P z m (cosh 0) =f— 7 E- / dip (cosh 9 + sinh 9 cos ip) 1 cos mip 

J o 


(l + m)! 


E 


i 


77 ^ fc!(Z — k)\ 


e (i-2k)8 / d^cos^-^^sin^^e 2 ^, 

















we find that 


/ (ji + m)P m _ i (cosh 9) P r _ 2 (cosh 9) \ / (j' + m)P n !_ ? (cosh#) (j' — m + 1)P™ ' ? (cosh#) 

3i 2 a Ji 2 J Ji+2 

\ — (ji — m + l)P m ! 2 (cosh 0) P'A 2 (cosh 9) ) \ —P™ + ? (cosh 9) P™~ l ~ 1 2 (cosh 9) 

\ 3 i~r 2 Ji+2 / \ h~ 2 Ji+2 

(j,+m)!(j' + m)! ^ 1 A 1 

^ to fc! (* + 3 - fc ) ! to fc'! t + I - fc') ! 


where 


x / dp dp' cos 2ji ~ 2k psin 2k pe 2imlp cos 2j '^- 2k ' p' sin 2k ' p'e 2imv ' 


x exp ( - H- m + — rii+i + qt — 2k - 2k! 9 

\ a 2a 2a 


A =--- 7 e"* (ji + k) U + fc') - e^+v'A , 

cosipcosp' \ 2 / \ 2 / V / 

B = (ji + i - fc) [(j' - m + 1) + (j[ + m + 1) e^+^1 , 

\ 2 J cos 99 J 

C = - (j( + \ - fc') ^ [(ji - m + 1) + (j, + to + 1) e^'>l , 

\ 2 J cos (p f J 

D = — e e cos p cos </?' (ji — to + 1) (j' — to + 1) e _l ^ +v ) — (j, ; + to + 1) (j' + to + 1) e^ ¥>+¥ ’ ^ . 


Making a change of variables 


‘P / <P 

^ = 71 ’ * = vr 


we have 


dip cos 2ji_2fc y> sin 2fc pe 2im,p 


l k+ i J - c 


Z fe+ 5 J -0 


/ \ 2 k 

dpp 2k f 1 - t) 6XP 


exp (-( 2 i + 2 ni- 2 ^( 1 - + ^ 


dpp 2k (1 + Pi, 2 ) exp -p 2 H-A- + A,i + A, 


Here and in the following, we denote by <¥,.1 and X.^i terms of order yfe and e respectively. Similarly, we have 


7 / 2i'—2k' / • 2k' / 2imip' 

dip cos <p sin ip e * 


l k ' + 5 J -c 


dp’p nk ' ( 1 - ) exp 


2irrup' 

~^/T 


2b, b 


exp rt + o (n * + n ' +l) + 2 ® - 2ft# J l^r + 12 P 


<£' 2 ^' 4 


J dp'p' 2k (1 + Pi, 2 ) exp 1' 2 + 2lT Jf + C h 1 + Ci, 2 ^ ■ 


On the other hand, making a change of variables 


where 


— 9 + 9q, 


sinh#o = 
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we have expansions of the form 


, ,3 /I — r 2ib 

B —Z 2 (1 + £^1 + £^ 2 ) j 
C ~ — -^-Z 2 (l+J-'j j i+J : 'j j 2)) 


D ~l 


2 / 1 + t 26^Z((^ + <p ; ) | m 


1 — t a 


v 7 ! 


Z>Z 


and 


nOO 

/ dO sinh 9e~ KL cosh 9 exp 
Jo 

- f°° dd sinh (e -+ 0„) e- (1+£)aJCOsh (^ +00 ) exp 


Z a +1 b . 

—I-- —Tii + —Tij+i + — 2A: — 2A; 

a 2a 2a 


l a +1 b . 

- H —-—rii + —n i+ 1 + qi -2k- 2k 
a 2a 2a 


(0 + do) 


1 + T 


2a ^ 4a ' 4a ' 


Vl -r 2 V 1 _ r y 

x (1 + ft,i + ft, 2 ) exp ('Hi, 1 + Hi, 2 ) • 


7/1 / ^ ^ /i 2 ^ 

dO exp--—0- 

ar 2ar ar 


a +1 6 . 

n, + -n i+1 + ft ) ^ 


Using Stirling’s asymptotic expansion for gamma functions, we also obtain an asymptotic expansion 

1 \!(ji - m)\ (, j[ - m)\(ji + to)! {j'i + ™)\ 




Collecting the terms, we have 


Vl. 


2 (b 


1 (b 


k'+i 


k' — h 


l k +«- 1 exj>l^+l i , 1 +l i ,2 ) . 


1 + T 


2S- + Tt-"i+4S: n i+l-l--5 L 


\/l - T 2 V 1 - r 

x J dp0 lk (1 + £?i, 2 ) exp ^ + A,i + A, 2 ^ 

x [ dp'p ,2k (1 + Vi 2 ) exp (--p’ 2 + + Ci p + C i)2 

J -00 V a VZ 


^ ^ fc!Zc'! 


/c—0 A;'—0 


/ J_ 

ar 2 ar 


si ( a - 1-1 

ar 


dOex p ( - 77—0 2 - — + ( -^~ n i + 7TT n *+i + Qi ) 0 


b 

2a 


x (1 + Gi, 1 + Gi, 2 ) exp {Ui'i + Ui, 2 ) exp ( — +Ij,i +X ii2 


m 


/ / 1 — t- 

1 +T sfl 


1 + £i, 1 + £j,2 


\ 


1 n , t- , t~ \ ,1 + r fi(<p + <f/) , m 

^-(1 + ^,1 + ^) 7i + u)) 

Now we expand the terms to obtain terms of order yfe and e respectively, and then perform the summation over k 
and k'. After summation, the integrations over <p and ip' are standard, followed by the integration over 6. We obtain 
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an expression of the form 


2 ar (1 + r 


7 tI \ 1 — r, 

X (1 + J7i,l + Hi, 2 ) 
Interchanging ni and ?Zj+i gives 


n ‘+TS”‘+l+T?-+3 

/ l 


exp - 


\ ar 

^2^ 

1 + 


21 \ aT"- + Ta n '+ l + * 


/ 

— (1 + /C^i + /C^) — \/l — T 2 


m 
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+ l m 


2ar /1+ r 


2a ' 4a "* + 1^4a"*^ 2^2 


irl \1 — T, 
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m 
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TO 
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Now for Tj+|f and rjgf, 
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,/j' + l(M 

1 — 2 —-- 

If (biz) If (biz) 

Kj'^bu) Kj'. + i(buj) 

Kj> ( biz ) 


we use Debye asymptotic expansions of modified Bessel functions 


Iv{vz) 
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ui (t(z)) 
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where 


We find that 
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. t 5 1 3 
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Ij' t (biz) 
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Kf(biz) 



TO 2 \ 

2w 7 y 
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Here 7\i is real and 7^,2 is complex. After expansion and simplification, we find that 


, /T+’ S . 0 


0 T+’ s „ 
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' — l- 
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PXP I - 7 - - - m + —g—K ? + <,) + —2i— ».".+> 
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For T^V’i and T> irj ^’ 2 , we have similar expansions 
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ifj 4 +i(aw) 


1 + i 


V(4( 


1 + Vi,! + Vi,2 


(«w) 


Finally, we obtain an expansion of the form 
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The imaginary part has been omitted since it is easy to verify that it would not contribute in the trace. 
Substituting into (O, we find that the term of order yV is zero and we have 


T?Cas — 
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2 S+1 7T 2 ta 


V (~l) a ' 

^ s +1 


dr 


s=0 


VT^ 


T z J 0 


d/ / 2 / dm / dn\ ... 


drix 


' — OO <7 —OO 


x exp 


2e(s + 1)1 (5 + l)m 2 br 


ar 

s —1 s 


blr 


~ - m+if 


2=0 


1 + ]T X] >v,.vvv, • X vv.2 • 

2=0J=2+l 


ii±il 


i—0 


Ibr J 


+ (s + 1)(^2 + 7 / 2 ) I ■ 


The integrations over Ui, in, l and r are standard and we obtain 


„ hcba 2 VV (-1 ) s+1 V / . l-(s + l) 2 \ \ 

(TTlFV r 3a6 ) + --j 

7 7r 3 Her at b 


1 - 


1 20 
3 _ Vr 2 


d d 


2880d 2 (r A + r B ) 


rA + r B 

















































One can easily verify that the leading order term 
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Ecus ~ 


7 hcr A r B 

28807T d 2 (r a + r B ) 


coincides with the proximity force approximation. Since 


1 

3 


20 

77T 2 


= 0.0438, 


the sign of the next-to-leading order term can be the same as or different from the leading term depending on the 
ratio of r a to r B - 

We can also compare this result to the results of scalar fields and electromagnetic fields [Hcj: 


Ecas 


Z^N 

^Cas 


Ecas 


■n^hcrAfB 
1440d 2 (r^i + tb) 
7 t 3 hcrArB 
1440d 2 (?y4 + rs) 
TT 3 hcrAl~B 

720 cP(r A +r B ) 


d 1 / d d 
1 T-t - (-1- 

r A + r B 3 \r a r B 


1 - 


1 - 


ta +r B 
d 


+ 


40 

— 2 


VA + Tb 


+ 


TV 

20 

7T 2 


d_ 

ta 

d 


d 

tb 

d 


ta r B 


For the case of two Dirichlet spheres, the sign of the next-to-leading order term can also be the same as or different 
from the leading term depending on the ratio of to r B - However, for the case of two Neumann spheres or two 
perfectly conducting spheres, the sign of the next-to-leading order term is always different from the leading term. 

It is also interesting to note that for scalar field, fermionic field as well as the electromagnetic field, the ratio of the 
next-to-leading order term to the leading order term always contain the term 

d 

rA + r B 


with same coefficient. It looks like this is a universal term that does not depend on boundary conditions. 


V. CONCLUSION 


In this work, we consider the Casimir interaction between two spheres that results from the vacuum fluctuations of 
a massless Dirac field with MIT-bag boundary conditions. The sphere-sphere Casimir interaction due to scalar fields 
and electromagnetic fields has been well-understood. However, to the best of our knowledge, we are the first one 
that investigates the fermionic interaction between two spheres. Using our formalism in 0 , we derive the functional 
representation of the Casimir interaction energy. The most technical part in the derivation is the computation of 
the translation matrix that relates the fermionic spherical waves in two different coordinate systems. We tackle the 
problem using the operator approach as in flfl It? . The result can be considered as an important byproduct of this 
work. 

From the formula of the Casimir interaction energy, we compute the large separation and small separation asymp¬ 
totic behaviors. As usual, the large separation asymptotic behavior is very easy to compute since it only involves a 
few terms that can be computed explicitly. We showed that when the separation between the spheres is large, the 
Casimir interaction energy behaves like 


Ecas 


Qhcr\r^ 
nL 5 ’ 


where L is the separation between the centers of the two spheres, and rA and r B are the radii of the spheres. We 
note that the order of the Casimir interaction energy is L -5 , which is intermediate between the order of the Casimir 
interaction energy between two Dirichlet spheres (L~ 3 ) and the order of the Casimir interaction energy between two 
Neumann spheres (L~ 7 ). 

In the small separation limit, we compute analytically the asymptotic behavior of the Casimir interaction energy 
up to the next-to-leading order term and find that 


Eca 


7 tt 3 hcrArB 

2880 d 2 (r A +r B ) 


1 - 


rA + r B 


1 20 
3 _ 77T 2 


d_ 

rA 


d 

r B 
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Here d is the distance between the two spheres. The leading order term agrees with the proximity force approximation. 
We also note that the term 

_ d 
ta + r B 

which appear in the ratio of the next-to-leading order term to the leading order term is ’universal’ among different 
types of quantum field. 

From the leading terms of the small separation and large separation asymptotic behaviors, we find that the Casimir 
interaction force is an attractive force in these two regimes. We believe that the force is always attractive at all 
separations, but it is not easy to deduce this from the expression for the Casimir interaction energy. 
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